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Pure Stationary States of Open Quantum Systems 



m 
O 
O 

> 
O 



> 



m ■ 
o ■ 

-i— > 



13 



Vasily E. Tarasov 

Skobeltsyn Institute of Nuclear Physics, Moscow State University, Moscow 119992, Russia^ 

Using Liouville space and superoperator formalism we consider pure stationary states of open 
and dissipative quantum systems. We discuss stationary states of open quantum systems, which 
coincide with stationary states of closed quantum systems. Open quantum systems with pure 
stationary states of linear oscillator are suggested. We consider stationary states for the Lindblad 
equation. We discuss bifurcations of pure stationary states for open quantum systems which are 
quantum analogs of classical dynamical bifurcations. 



I. INTRODUCTION 

The open quantum systems are of strong theoretical 
interest. As a rule, any microscopic system is always em- 
bedded in some (macroscopic) environment and therefore 
it is never really closed. Frequently, the relevant environ- 
ment is in principle unobservable or it is unknown 0] . 
This would render the theory of open quantum systems a 
fundamental generalization of quantum mechanics 0, 0] . 

Classical open and dissipative systems can have regu- 
lar or strange attractors 0, 13 • Regular attractors can be 
considered as a set of (stationary) states for closed clas- 
sical systems correspondent to open systems. Quantiza- 
tion of evolution equations in phase space for dissipative 
and open classical systems was suggested in Refs. 0H|. 
This quantization procedure allows one to derive quan- 
tum analogs of open classical systems with regular at- 
tractors such as nonlinear oscillator ■ In the papers 
0, were derived quantum analogs of dissipative sys- 
tems with strange attractors such as Lorenz-like system, 
Rossler and Newton-Leipnik systems. It is interesting to 
consider quantum analogs for regular and strange attrac- 
tors. The regular "quantum" attractors can be consid- 
ered as stationary states of open quantum systems. The 
existence of stationary states for open quantum systems 
is an interesting fact. The condition given by Davies 
in Ref. [T(i| defines the stationary state of open quan- 
tum system. An example, where the stationary state is 
unique and approached by all states for long times is con- 
sidered by Lindblad pd| for Brownian motion of quan- 
tum harmonic oscillator. In Refs. 0EI3 Spohn 
derives sufficient condition for the existence of a unique 
stationary state for the open quantum system described 
by Lindblad equation [l^ . The stationary solution of the 
Wigner function evolution equation for an open quantum 
system was discussed in Refs. 0,0. Quantum effects 
in the steady states of the dissipative map are considered 
in Ref. [H|. 

In this paper we consider stationary pure states of some 



open quantum systems. These open systems look like 
closed quantum systems in the pure stationary states. 
We consider the quantum analog of dynamical bifurca- 
tions considered by Thompson and Lunn 0] for classi- 
cal dynamical systems. In order to describe these sys- 
tems, we consider Liouville-von Neumann equation for 
density matrix evolution such that this Liouville gener- 
ator of the equation is a function of some Hamiltonian 
operator. Open quantum systems with pure stationary 
states of linear harmonic oscillator are suggested. We 
derive stationary states for quantum Markovian master 
equation usually called the Lindblad equation. The sug- 
gested approach allows one to use theory of bifurcations 
for a wide class of quantum open systems. We consider 
the example of bifurcation of pure stationary states for 
open quantum systems. 

In Section II, we introduce pure stationary states for 
closed quantum systems and some mathematical back- 
ground is considered. In the Section III, we consider 
Liouville-von Neumann equation for an open quantum 
system and pure stationary states for this equation. In 
Section IV, simple examples stationary states for open 
quantum systems are considered. In Section V, we study 
some properties of quantum system to have dynamical 
bifurcations and catastrophes. In the Section VI, we 
suggest an example of quantum system with fold catas- 
trophe. Finally, a short conclusion is given in Section 
VII. In Appendix, the mathematical backgroung (Liou- 
ville space, superoperators) is considered. 



II. PURE STATIONARY STATE 

In the general case, the time evolution of the quantum 
state \pt) can be described by the Liouville-von Neumann 
equation 



dt 



\pt)=Mfh), 



(1) 
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where A is a Liouville superoperator on Liouville space, 
\p) is a density matrix operator as an clement of Liouville 
space. For the concept of Liouville space and superop- 
erators see the Appendix and Refs. [2(j-[37|]. For closed 



2 



systems, Liouville superoperator has the form 



A = --{L H -R H ) or A = L H , 



(2) 



where H = H(q,p) is a Hamilton operator. If the Li- 
ouville superoperator A cannot be represented in the 
form Pjl. then quantum system is called open, non- 
Hamiltonian or dissipative quantum system [8ft IstI l38j ] . 
The stationary state is defined by the following condition: 



Mpt) = o. 



(3) 



For closed quantum systems (J2Jl, this condition has the 
simple form: 



L H \pt) = Rh\pi) or L H \p t ) = 



(4) 



In the general case, we can consider the Liouville super- 
operator as a superoperator function 7, 8, 36] : 

A = A(L X ,L+) or A = A(L X ,R X ), 

where X is a set of linear operators. For example, X — 
{q,p,H} or X — {Hi, ..,H S }. In the paper we use the 
special form of the superoperator A such that 

s 

A = - Rh) + ^F k N k (L H ,R H ), 

k=l 

where N k (Ln, Rh) are some superoperator functions and 
F k is an arbitrary nonzero superoperator. 

It is known that a pure state \pv) is a stationary state 
of a closed quantum system Q, (J2J, if the state \py) is 
an eigenvector of the Liouville space for superoperators 
Lh and Rh- 



Lh\p<$) = \p*)E, Rh\p*) = \p*)E. 



(5) 



Equivalently, the state \pqr) is an eigenvector of superop- 
erators L H and L H such that 

l h\p*) = \p<sr)E, L H \p-q,) = \p 9 ) ■ = 0. 

The energy variable E can be defined by 

E = (I\L H \pv) = {I\Rh\pv) = (I\L H \p*). 

The superoperators Lh and Rh for linear harmonic 
oscillator are 

f 1*2 rnu) 2 f 2 - 1 « 2 mu> 2 - 2 

L H = —L -\ — L R H = —R p -\ — R q - (6) 

2m 1 2 H 2m 1 2 H 

It is known that pure stationary states p^ n = p\ of 
linear harmonic oscillator JJ]J exists if the variable E is 
equal to 



III. PURE STATIONARY STATES OF OPEN 
SYSTEMS 

Let us consider the Liouville-von Neumann equation 
(|T|) for the open quantum system defined of the form 



dt 



) = -j(L H - R H )\pt) + Y,F k N k {L H ,R H )\p t ) 



k=l 



(8) 

Here F k is some superoperator and N k (Ln , Rh), where 
k = 1, s, are superoperator functions. 

Let \p<Sf) is a pure stationary state of the closed quan- 
tum system defined by Hamilton operator H . If Eqs. (J5J 
are satisfied, then the state \py) is a stationary state of 
the closed system associated with the open system |JSJ 
and is defined by 



j t \Pt) = -^H-R H )\p t ). 



(9) 



If the vector \p^) is an eigenvector of operators Lh and 
Rh, then the Liouville-von Neumann equation (JHJ for the 
pure state \p-q,) has the form 



dt 



IP* 



)=E EM N k (E,E), 



k=l 



where the function N k (E, E) are defined by 

N k (E,E) = {I\N k {L H ,R H )\p*). 

If all functions N k (E, E) are equal to zero 

N k {E,E) = 0, (10) 

then the stationary state \p^) of the closed quantum sys- 
tem (0 is the stationary state of the open quantum sys- 
tem ©. 

Note that functions N k (E, E) are eigenvalues and \py) 
is the eigenvector of superoperators N k {LH,R.H), since 

N k (L H ,R H )\p*) = \p*)N k {E,E). 

Therefore stationary states of the open quantum sys- 
tem (jHJ are defined by zero eigenvalues of superoperators 
N k (L H ,R H ). 



IV. OPEN SYSTEMS WITH OSCILLATOR 
STATIONARY STATES 

In this section simple examples of open quantum sys- 
tems (jSJ are considered. 

1) Let us consider the nonlinear oscillator with friction 
defined by the equation 



E n = -Tvuj{2n+l) 



j t Pt = -J^HiPt] - 7^rP[q 2 ,P 2 Pt + Ptp 2 }, (11) 
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where the operator H is the Hamilton operator of the 
nonlinear oscillator: 



H = 



p mil q 7<7 



2m 2 2 
Equation can be rewritten in the form 
d 



dt 



\Pt) = L H \p t )+ 



7 A 



wW-w^M* (12) 



where A = f2 2 — oj 2 and the superoperator is defined 
for the Hamilton operator H of the linear harmonic os- 
cillator by Eqs. J5J) and ©. Equation (jT2"|) has the form 
©, with 

N(L H ,R H ) = ^{L H +R H )-^L I , F = 2mpL- 2 . 
In this case the function N(E, E) has the form 



N(E, E) = E 



2/3' 



Let 7 = (5m 2 u 2 . The open quantum system ifTTl) has one 
stationary state of the linear harmonic oscillator with en- 
ergy E n = (huj/2)(2n+l), if A = 2j3huj(2n+l), where n 
is an integer non-negative number. This stationary state 
is one of the stationary states of the linear harmonic os- 
cillator with the mass m and frequency uj. In this case 
we can have the quantum analog of dynamical Hopf bi- 
furcation pjjl3^|. 

2) Let us consider the open quantum system described 
by the time evolution equation 

j t \pt) = L H \ Pt ) + L-cos(^L+) \p t ), (13) 

where the superoperator is defined by formulas (|2")l 
and ©■ Equation (|13fl has the form (JSJ if the superop- 
erators F and N(Lu, Rh) are defined by 

F = --r{L q - R q ), 



N(L H ,R H ) 



cos(^-(Lh + Rh) ) = 



E 

m=0 



2 m 



(14) 



The function N(E, E) has the form 

w( ,, E) = ™(f) = £^( £o 



The stationary state condition l|10[l has the solution 

£=|(2n + l), 

where n is an integer number. If parameter Sq is equal 
to tiu), then quantum system (|13|l and l|14|) has pure sta- 
tionary states of the linear harmonic oscillator with the 
energy Q. As the result, stationary states of the open 
quantum system (| 1 31) coincide with pure stationary states 
of linear harmonic oscillator. If the parameter £o is equal 
to hcu(2m + 1), then quantum system l(T3"|l and have 
stationary states of the linear harmonic oscillator with 
n(k, m) — 2km + k + m and 

E n{k , m ) = ^(2fc + l)(2m+l). 

3) Let us consider the superoperator function 
Nk{Ln, Rh) in the form 
1 



N k {L H ,R H ) 



2% ^ 



11 T1* (IT 71 D m T n + 

v kn v km\ z - u H It H ^ H 



-R^ +m ) 



and all superoperators F k are equal to Lj. In this case, 
the Liouville-von Neumann equation W ca n be repre- 
sented by the Lindblad equation [3|| l40l l4l| : 

d 



dt 



\pt) 



-{L H - R H )\p t )- 



in 



2L Vk R v , - L V L V , - R v ,R Vk \ Pt ). (15) 



(16) 



with linear operators V k defined by 



If \p<st) is a pure stationary state (JHJ, then all functions 
N k (E,E) are equal to zero and this state p*) is a sta- 
tionary state of the open quantum system (|15|l . 
If the Hamilton operator H is defined by 

™, ,2 



H 



1 



A 



-V 



(qp+pq), 



(17) 



2m 2 2 

then we have some generalization of the quantum model 
for the Brownian motion of a harmonic oscillator consid- 
ered in [llj . Note that in the model [ll| operators V k are 
linear V k = a k p + b k q, but in our generalization i|15|) and 
(|16|l these operators are nonlinear. For example, we can 
use 

V k = a k p + b k q + c k p 2 + d k q 2 + e k (qp + pq). 

The case c k — d k = e k = is considered in Ref. [Tl| . 
Let real parameters a and f3 exist and 



b k = aa k , c k = Pa k , d k 



m 2 uj 2 (3a k , 



e k = m\(3a k . 



In this case, the pure stationary states of the linear os- 
cillator (fP7)l exist if ui > A and the variable E is equal 
to 



E n = hiuj{2n + 1)^/1 - X 2 /uj 2 . 
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V. DYNAMICAL BIFURCATIONS AND 
CATASTROPHES 

Let us consider a special case of open quantum systems 
(JSJl such that the vector function 

N k (E,E) = (I\N k (L H ,R H )\p), 

be a potential function and the Hamilton operator H can 
be represented in the form 



N 



k=l 



In this case we have a function V(E) called potential, 
such that the following conditions are satisfied: 

dV(E) 



dE k 



N k (E,E). 



where E k = (I\L Hk \p) = (I\R Hk \p). If potential V(E) 
exists, then the stationary state condition (|10fl for the 
open quantum system ((SJ) is defined by critical points 
of the potential V(E). If the system has one variable E, 
then the function N(E, E) is always a potential function. 
In general, the vector function N k (E, E) is potential, if 

dN k (E,E) dNi(E,E) 



dE x 



dE k 



Stationary states of the open quantum system (JSJ) with 
potential vector function N k (E,E) is depend by critical 
points of the potential V(E). It allows one to use the 
theory of bifurcations and catastrophes for the paramet- 
ric set of functions V(E). Note that a bifurcation in a 
vector space of variables E = {E k \k = 1, s} is a bifur- 
cation in the vector space of eigenvalues of the Hamilton 
operator H k . 

For the polynomial superoperator function 
N k (L H ,R H ) we have 

N n 

N k {L H ,R H ) = £ £ ai%L^R n H m . 

n—0 m— 

In general case, m and n are multi-indices. The function 
N k (E, E) is a polynomial 

N 



N k {E,E) = Y. a n )E7] 



n=0 



where the coefficients affi are defined by 



»(*) 



V a {k) 

m=0 



We can define the variables xi = Ei — a\ (I = 1, s), 
such that functions N k (E, E) = N k (x + a,x + a) have no 
the terms x? -1 . 



N 

N k {x + a,x + a) = J2 a( n } (x + a (k) ) n 

n=0 



EE< 

n— m=0 



„0). 



m\(n — my. 



x m {a {k) ) n -' m . 



If the coefficient of the term x? 1 1 is equal to zero 

(fc) _^!_ ( fe ) ( k) (fe) ( k ) (fe) Q 

1 (rii — 1)! "ni-l — "rti U ' U Z u ' 

then we have the following coefficients 



(fe) 

(fe) _ a m -l 
a l ~ (fe) ■ 

TllOCfn 

If we change parameters , then an open quantum 
system can have pure stationary states of the system. For 
example, the bifurcation with the birth of linear oscillator 
pure stationary state is a quantum analog of dynamical 
Hopf bifurcation 0, |3^| for classical dynamical system. 

Let a vector space of energy variables E be a one- 
dimensional space. If the function N(E, E) is equal to 

n-l 

N(E, E) = ±a n E n + (*jE j n > 2, 
2=1 

then the potential V(x) is defined by the following equa- 
tion: 

n— 1 

V(x) = ±x n+1 + aj x j n > 2, 
2=1 

and we have catastrophe of type A± n . 

If we have s variables Ei, where I = l,2,...,s, then 
quantum analogs of elementary catastrophes A± n , D± n , 
E±6, Ej and Eg, can be realized for open quantum sys- 
tems. Let us write the full list of potentials V(x), which 
leads to elementary catastrophes (zero modal) defined by 
V(x) — Vq(x) + Q(x), where 

n-l 

A± n : V (x) = ±xl +1 + a 3 x{ n > 2, 

2=1 

n — 3 n-l 

D± n : V {x)=xix 2 ±x r >- 1 + Y,a J xi+ £ x{- ( " " 



j=l j=n-2 
2 5 



E±e ■ Vq(x) = (x\ ± x\) + a,jX2 + cijXi: 



2=1 2=3 
4 6 



E? : Vq(x) — x\ + X1X2 + /J ajX J 2 + QjXiC 



.2-3 



2-5 
2 i 



2=1 2=5 
3 7 

Eg : Vq(x) = x\ + xl + ^ djxi + ^ ajXixP 2 " A . 

2=1 2=4 

Here Q(x) is the nondegenerate quadratic form with vari- 
ables X2, X3, x s for A± n and parameters £3, x s for 
other cases. 
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VI. FOLD CATASTROPHE 

In this section, we suggest an example of the open 
quantum system with catastrophe A 2 called fold. 

Let us consider the Liouville-von Neumann equation 
(|5|) for a nonlinear quantum oscillator with friction, where 
multiplication superoperators Lh and Rh are defined by 
Eq. © and superoperators F and N(Lh, Rh) are given 
by the following equations: 



F 



-2L-L+ 



(18) 



N(L H ,R H ) =a Q L+ + ai L++a 2 (L+) 2 , (19) 

In this case, the function N(E, E) is equal to 

N(E, E) = a + a x E + a 2 E 2 . 

A pure stationary state ) of the linear harmonic os- 
cillator is a stationary state of the open quantum system 
(^HJ), if N{E, E) = 0. Let us define the new real variable 
x and parameter A by the following equation 



E 



2a 2 



A 



4aoa 2 — a\ 



4a! 



Then we have the stationary condition N(E, E) = in 
the form x 2 — A = 0. If A < 0, then the open quan- 
tum system has no stationary states. If A > 0, then we 
have pure stationary states for a discrete set of parame- 
ter values A. If the parameters ai,a 2 and A satisfy the 
following conditions 

ai . , 1 , m. . _ 2 ? m2 

= huj(n-\ 1 ), A = h uj — , 

2a 2 1 2 2 h 4 ' 

where n and m are non-negative integer numbers, then 
the open quantum system H18|) and 1)19(1 has two pure 
stationary state of the linear harmonic oscillator. The 
energies of these states are equal to 

E n = TiLo(n + -) , E n+m = %uj(n + m+ -). 



VII. CONCLUSION 

Open quantum systems can have pure stationary 
states. Stationary states of open quantum systems can 
coincide with pure stationary states of closed (Hamilto- 
nian) systems. As an example, we suggest open quantum 
systems with pure stationary states of linear oscillator. 
Note that using Eq. J&J, it is easy to get open (dissipa- 
tive) quantum systems with stationary states of hydro- 
gen atom. For a special case of open systems, we can use 
usual bifurcation and catastrophe theory. It is easy to de- 
rive quantum analogs of classical dynamical bifurcations. 
Physical instances with bifurcation behavior can be real- 
ized in quantum optics |42| and deep inelastic collisions 

m. 



Open quantum systems with two stationary states can 
be considered as qubit. It allows one to consider open 
n-qubit quantum system described by Eq. Q as a quan- 
tum computer with pure states. In general, we can con- 
sider open n-qubit quantum systems as quantum com- 
puter with mixed states 0,^3- A mixed state (operator 
of density matrix) of n two-level quantum systems (open 
or closed n-qubit system) is an element of 4 n -dimensional 
operator Hilbert space (Liouville space). It allows one 
to use quantum computer model with four-valued logic 
|44l |4(| . The quantum gates of this model are real, com- 
pletely positive, trace-preserving superoperators that act 
on mixed state. Bifurcations of pure quantum states can 
used for quantum gate control. 
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Appendix 

For the concept of Liouville space and superoperators 
seeRefs. IH-Hl 



Liouville space 

The space of linear operators acting on a Hilbert space 
H is a complex linear space TC. We denote an element A 
of TL by a ket-vector \A). The inner product of two ele- 
ments \A) and \B) of H is defined as (A\B) = Tr(A^B). 
The norm = y/(A\A) is the Hilbert-Schmidt norm 
of operator A. A new Hilbert space TL with the inner 
product is called Liouville space attached to Tt or the 
associated Hilbert space, or Hilbert-Schmidt space po|- 

M 

Let {\x >} be an orthonormal basis of Tt: 



< x\x' >— 6(x — x ) , 



dx\x >< x\ = I. 



Then \x,x') = \ \x >< x'\) is an orthonormal basis of the 
Liouville space Tt: 



{x,x'\y,y') = 5{x - x')8{y - y'), 



dx J dx'\x, x')(x, x'\ = I. 
For an arbitrary element \A) of Tt we have 
\A)= J dx J dx'\x,x')(x,x'\A) 

where 

(x,x'\A)=Tr((\x ><:z'|)U) = 



(20) 



(21) 
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= Tr(\x' >< x\A) =< x\A\x' >= A(x,x'), 

is a kernel of the operator A. An operator p of density 
matrix (Trp = 1, p' = p, p > 0) can be considered as an 
element |p) of the Liouville space Ti. Using (|21[) . we get 



Left superoperators are defined as Lie and Jordan 
multiplication by 



1 



1 



L A B = —(AB - BA), L+B = -(AB + BA). 



\p) = J dx J dx'\x,x')(x,x'\p) , 
where the trace is represented by 

(I\p) = Trp = dx (x, x\p) = 1. 



(22) 



Superoperators 

Operators that act on Ti are called superoperators and 
we denote them, in general, by the hat. 

For an arbitrary superoperator A on W, which is de- 
fined by 



A|A) = \k(A)), 



we have 



The left superoperators L A and right superoperators R A 
are connected by 



L A - -R A , L\ - R\. 



An algebra of the superoperators L A is defined Q by the 
following relations 
1) Lie relations 



l ab — l a l b L B L A . 



2) Jordan relations 



f + -t- T+T + T+ 4- T + T + T + 

L (AoB)oC + ^B^C^A + ^A^C^B 



T+ t+ 4- T + f+4-f + f + 
^AoB^C + ^BaC^A + ^AoC^B' 



(x,x'\k\A) = dy dy'(x,x'\A\y,y')(y,y'\A) = 



L (AoB)oC + L B L C L A + L A L C L 1 



dy / dy'A(x,x',y,y')A(y,y'), 



T + T+ 4- T + T+ 4- T + T + 
^C^AoB + ^B^AoC + L A L BoC 



where A(x, x', y, y') is a kernel of the superoperator A. 

Let A be a linear operator in the Hubert space Ti. We 
can define the multiplication superoperators La and Ra 
by the following equations 

L A \B) = \AB) , R A \B) = \BA). 

The superoperator kernels can be easy derived. For 
example, in the basis \x,x') we have 

(x,x'\L A \B) = [dy [ dy'(x,x'\L A \y,y')(y,y'\B) = 



T + T+ J- T + T+ 4- T + T + 

^C^AoB + ^B^AoC + ^A^BoC 



T + T+ 4- T + T+ 4- T + 



3) Mixed relations 
L 



AB — A B ^B^Ai 



L AoB - L A L B + L B L A' 



dy I dy'L A (x 1 x' ,y,y')B(y,y'). 



Using 



(x,x'\AB) =< x\AB\x' >= 

= J dy J dy' < x\A\y >< y\B\y' >< y'\x' >, 

we get kernel of the left multiplication superoperator 
(L A )(x,x',y,y') =< x\A\y >< x'\y' >= A(x,y)5(x' -y'). 



Laob — L A L B —L b L a , 



T + T+ f + f+ — f- 
^B^A _ ^A^B — ~^ A Bi 



where 



A- B = ^-(AB-BA), AoB = -(AB + BA). 
in 2 
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